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Abstract. Linear cosmological perturbations of a large class of modihed gravity and dark 
energy models can be unihed in the effective held theory of cosmic acceleration, encompassing 
Horndeski scalar-tensor theories and beyond. The fully available model space inherent to this 
formalism cannot be constrained by measurements in the quasistatic small-scale regime alone. 
To facilitate the analysis of modihcations from the concordance model beyond this limit, we 
introduce a semi-dynamical treatment extrapolated from the evolution of perturbations at a 
pivot scale of choice. At small scales, and for Horndeski theories, the resulting modihcations 
recover a quasistatic approximation but account for corrections to it near the Hubble scale. 
For models beyond Horndeski gravity, we hnd that the velocity held and time derivative of 
the spatial metric potential can generally not be neglected, even in the small-scale limit. We 
test the semi-dynamical approximation against the linear perturbations of a range of dark 
energy and modihed gravity models, hnding good agreement between the two. 
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1 Introduction 

Identifying the nature of the late-time accelerated expansion of our Universe is a prime en¬ 
deavour of cosmology. Instead of a vacuum energy in the form of an inexplicably small 
cosmological constant, cosmic acceleration may be driven by the contribution of a dark en¬ 
ergy or a modification of gravity. Conceptually, gravitational physics is not understood in the 
ultraviolet and a more fundamental theory of gravity may give rise to a remnant in the infrared 
that could cause the effect. The simplest potential contribution is a single low-energy effective 
scalar degree of freedom. A plethora of modified gravity models have been proposed based 
on the prospects of such a field [1-3] and a more systematic approach to explore their cosmo¬ 
logical implications has become a necessity. In the past few years, extensive efforts have been 
conducted to develop a generalised formalism for the resulting modifications in the formation 
of structure [4-23]. The effective field theory (EFT) of cosmic acceleration [14, 18, 19, 21-28], 
or unified dark energy, provides such a framework and describes the evolution of the spatially 
homogeneous and isotropic background of our Universe and the perturbations around it for 
a large class of modified gravity and dark energy scenarios. 

Solving the coupled system of linearly perturbed modified Einstein equations for gen¬ 
eral dark energy and modified gravity models can be a time-consuming numerical challenge 
and may only offer little insight into the generic behaviour of the wide range of possible 
modifications. Contrary to the A Cold Dark Matter (ACDM) concordance model, where the 
evolution of the gravitational potential is scale independent, for modified models the field and 
conservation equations need to be re-evaluated at each scale of interest. To circumvent the 
consequent increase in computational demand, and to enable analytic calculations, frequently 
a quasistatic approximation is employed, which neglects time derivatives compared to spatial 
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derivatives and usually emphasises density perturbations compared to large-scale velocity cur¬ 
rents in the held equations. This can provide direct information about the relations between 
the huctuations without requiring the integration of the held and conservation equations. 
However, the contributions of time derivatives of the metric potentials and the velocities can 
become important near the Hubble scale, or near the sound horizon for models with sublu¬ 
minal sound speed in the dark energy component. These effects can limit tests of gravity 
and dark energy employing the quasistatic approximation and including observations at very 
large scales [20, 29, 30]. Importantly, one of the free time-dependent functions characterising 
the EFT model space does not enter the quasistatic computations [28] and, hence, cannot be 
constrained by measurements in the quasistatic regime. An analysis beyond the quasistatic 
limit is also required to lift degeneracies between the standard concordance cosmology and 
other subsets of the EFT model space [28]. Moreover, the cosmological background and the 
linear leading-order modihcations at small scales, characterised by a deviation in the Poisson 
equation and a gravitational slip between the metric potentials, only allow a measurement of 
three independent functions of time. In principle, a measurement of the linear propagation 
speed and friction of gravitational waves could be used to further restrict the available model 
space [31-34] but the required observations associated with a late-time cosmological modih- 
cation are currently lacking. However, since the EFT formalism contains more freedom than 
this, in order to fully constrain the available model space, next-to-leading-order contributions 
to these modihcations need to be measured. This implies the consideration of the velocity 
held and the time derivatives of the huctuations. Hence, a more general, yet efficient ap¬ 
proach beyond the quasistatic approximation that also allows for direct analytical operations 
to gain insights on the rich phenomenology of dark energy and modihed gravity models is 
highly desirable. 

In this paper, we adopt the unihed dark energy formalism of Ref. [22, 23], encompassing 
Horndeski scalar-tensor theories [35, 36] and models beyond it [37], and explore the linear 
cosmological perturbations that it describes. We extend previous work by deriving the energy- 
momentum conservation equations and a set of combined modihed Einstein equations where 
the scalar held huctuations are eliminated and which incorporate beyond-Horndeski models. 
We then introduce a semi-dynamical treatment extrapolated from the evolution of the pertur¬ 
bations at a chosen scale, the pivot scale, to facilitate the analysis of modihcations from the 
concordance model and its application to observational tests of gravity and dark energy. We 
describe the resulting effective deviation in the Poisson equation, the gravitational slip, and 
an approximated growth rate of matter density huctuations. This also includes a discussion 
on the choice of pivot scale and its impact on the computation of these relations. Finally, 
we numerically test the semi-dynamical approximation against the exact linear perturbations 
for a range of dark energy and modihed gravity models that are embedded in the Horndeski 
and beyond-Horndeski action. These include models with a kinetic contribution of the scalar 
held, a running Planck mass, kinetic braiding, deviations between the speed of gravitational 
waves and the speed of light, as well as beyond-Horndeski terms. We also compare our results 
against a quasistatic approximation. 

We begin by briehy reviewing the unihed dark energy formalism in Sec. 2. In Sec. 3, 
we then introduce the semi-dynamical approximation for linear cosmological perturbations of 
unihed dark energy, determining the time and scale dependence of the effective modihcations. 
We test the semi-dynamical approximation against numerical solutions of the exact modihed 
Einstein and energy-momentum conservation equations in Sec. 4 and discuss our results in 
Sec. 5. The modihed linearly perturbed Einstein equations and coefficients of the effective 
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modifications are given in the appendices for reference. 


2 Unified dark energy 


The effective field theory of cosmic acceleration [18, 19] provides a unified framework for de¬ 
scribing the expansion of a Friedmann-Lemaitre-Robertson-Walker (FLRW) background and 
perturbations around it in a generalised theory of gravity constructed from supplementing 
the tensor field with a single scalar field and with the matter sector obeying the weak equiv¬ 
alence principle. Adopting the unitary gauge, where the time coordinate is chosen in order 
to absorb the scalar field perturbation in the metric the unified treatment reduces to 
finding a generalised action composed of geometric operators that are invariant under time- 
dependent spatial diffeomorphisms with free time-dependent coefficients. The EFT action, 
up to quadratic order, becomes [18, 19] 


5 ^ y [ mR - 2A(t) - T{t)5g^^^ + Ml{t){5g^y - Mt{t)5g^HK^^ 


T<Sni [V^mj g^h 


( 2 . 1 ) 


where = SvrG with bare gravitational constant G and the speed of light in vacuum is 
set to unity. Furthermore, R and R^^^ denote the four-dimensional and spatial Ricci scalar, 
respectively, is the extrinsic curvature tensor, and describes the normal to surfaces 
of constant time, whereby 5 denotes perturbations with respect to the background. We have 
adopted the notation of Ref. [28] for the EFT coefficients. Moreover, we only considered the 
leading order contributions with respect to their mass dimension in Eq. (2.1) and neglected 
terms that introduce higher-order time derivatives in the Euler-Lagrange equations. 

Given the quadratic order, Eq. (2.1) fully describes the background evolution and lin¬ 
ear perturbations of an extensive class of scalar-tensor theories. In particular it encom¬ 
passes Horndeski theory [35, 36], which describes the most general local, Lorentz-covariant, 
and four-dimensional scalar-tensor gravity where the Euler-Lagrange equations are at most 
second-order in the derivatives of the scalar and tensor fields. In Horndeski gravity, all coef¬ 
ficients, except for m 2 , are used in Eq. (2.1) with the restriction that M| = — M| = 2M^. 
Healthy theories beyond the Horndeski action [37], introducing third-order derivatives in the 
modified Einstein equations but with a constraint equation ensuring a second-order equa¬ 
tion for the propagating scalar degree of freedom, are also embedded in Eq. (2.1), allowing 
M| 7 ^ 2M^. Finally, a nonzero m 2 is, for instance, introduced in Hofava-Lifshitz [38] gravity 
with the violation of Lorentz covariance. The concordance model phenomenology is recovered 
when H = 1 , A is a constant, and the remaining coefficients vanish, but it can be degener¬ 
ate with other choices of the EFT functions on quasistatic scales [28]. In general, the EFT 
action, Eq. (2.1), introduces nine time-dependent coefficients to which the statistically spa¬ 
tially homogeneous and isotropic background metric adds the scale factor a(t), or the Hubble 
parameter H{t). The Friedmann equations introduce two constraints such that the formal¬ 
ism is composed of eight free time-dependent functions, given the matter content and spatial 
curvature fco [18, 19, 28]. Importantly, the EFT function M^Kt) is not contributing in the 
quasistatic limit [28]. 

In the following, we study the linear scalar perturbations around the FLRW background 
metric in unified dark energy, restricting to a spatially flat (fco = 0 ) universe, otherwise 
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containing only pressnreless dnst Pm = 0. We fnrthermore restrict to models where = 
— M| and m 2 = 0 in the EFT action, Eq. (2.1). This encompasses Horndeski gravity, where 
additionally 2M^ = M 21 beyond-Horndeski models with the corresponding effective 
theory described by five and six free fnnctions of time, respectively. We adopt the Newtonian 
gange with the FLRW line element described by 

ds2 = _(i + 2T)dt2 + a^{t){l + 2$)dx2 (2.2) 


and matter ffnctnations are considered in the total matter gange with velocities and density 
finctnations Am aronnd the backgronnd matter density pm- However, instead of working 
with the effective field theory coefficients of Eq. (2.1), we choose to adopt a description 
that characterises the free time-dependent fnnctions of the nnified dark energy with more 
direct physical, and observational, implications. Snch a framework was developed in Ref. [22] 
describing the backgronnd and linear pertnrbation theory eqnivalent to Horndeski gravity. 
We adopt the slightly different notation of Ref. [23], thereby also employing their extension 
that incorporates beyond-Horndeski models. The formalism separates ont the expansion 
history H as the free fnnction determining the cosmological backgronnd, which relates to the 
EFT fnnctions H, T, and A via the Friedmann eqnations. Linear pertnrbations aronnd the 
backgronnd are then characterised by five additional free fnnctions of time a*, each implying 
a different physical effect in the properties of the nnified dark energy. 

• Kineticity a-^: The contribntion of a kinetic energy of the scalar field can give rise to 
a clnstering of the dark energy component at very large scales. It relates to Eq. (2.1) as 

T -H , 


• Planck mass evolution rate The evolution of the gravitational coupling gives 
rise to a gravitational slip between the metric potentials T and <h. This attributes an 
effective anisotropic stress to the dark energy component. In terms of the EFT functions 
of Eq. (2.1), the parameter is expressed as 


+ (M|)' 


(2.4) 


where primes denote derivatives with respect to In a here and throughout the paper. 
The Planck mass is + M|) with om = dlnM^/dlna and relates to the 

background expansion via H = a^dlnM"^/dt. 


• Braiding a-Qi The interaction of the scalar field and the metric through braiding, or 
mixing, of the kinetic contributions of these fields causes the dark energy component to 
cluster at small scales. The relation to the EFT framework is given by 


HQ.' + Mf 

2H{Q + M^y 


(2.5) 


• Tensor speed alteration aT. A deviation between the speed of gravitational waves 
and the speed of light contributes as an effective anisotropic stress and clustering of the 
dark energy component. The parameter relates to Eq. (2.1) via 


M| 

= -QTwr 


( 2 . 6 ) 
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Beyond-Horndeski term a-^: The introduction of third-order derivatives in the Ein¬ 
stein equations, however, with a second-order equation for the propagating scalar degree 
of freedom, lies outside of the model space of Horndeski gravity [37]. As will be discussed 
in Sec. 3, this causes velocity helds and the rate of evolution of the metric potentials to 
contribute to the anisotropic stress and clustering of the dark energy on small scales at 
leading order. The relation to the EFT formalism is provided by 


OH = 


2M2 - M| 


(2.7) 


Note that the functions aj can also directly be dehned from the quadratic expansion of the 
action composed in terms of geometric quantities invariant under spatial diffeomorphisms in 
the 3 -|- 1 Arnowitt-Deser-Misner (ADM) spacetime decomposition with uniform scalar held 
hypersurfaces of constant time [23]. 

In order to describe the cosmological perturbations implied by Eqs. (2.1) and (2.2), the 
time diffeomorphism t —>■ t -|- 7r(t, x) with huctuation of the scalar degree of freedom vr is 
applied to the action to restore its full four-dimensional covariance. We work in Fourier space 
and simplify notation by referring to perturbative quantities through their Fourier amplitudes 
of the plane waves with comoving wavenumber k, where in linear theory the phases factor out 
in the held equations. We refer to Ref. [23] or Appendix A for the resulting perturbed modihed 
Einstein equations. We follow Ref. [22] to combine the held equations in two modihed Einstein 
equations where contributions of tt, its derivatives, and T' are eliminated. This implies using 
vr" determined by the scalar held equation, Eq. (A. 5), in the time-space component Eq. (A.2) 
to get vr', which then is used in the traceless space-space component Eq. (A.3) to obtain vr. 
The time-time component Eq. (A.l) then determines T, while T' can be obtained from the 
time derivative of the traceless space-space component. Finally, we combine these results in 
the trace of the space-space component Eq. (A.4) to get the hrst combined modihed Einstein 
equation. 


+ Ai4> -|- A2k^^ — Pm ^AsAm + ^ ) 

(2.8) 

where the coefficients are given by 


^ _ H' Pu + PsiaB - an^kj^ 

H + (o-g _ an)'^k‘jj- 

(2.9) 

^ + /3i 5 + (as - anf‘clk‘]j 

^1 + (as - anykjj 

(2.10) 

1 ^16 + 

^ - 2772M2^, + («g_«g)2fc2^’ 

(2.11) 

^ 1 PiQ8 +PidioB - OiH)‘^kjj + auenkjj 

2772M2 + («g _ 

(2.12) 

with kn = k/(aH) and 


i^i = l^i + 

(2.13) 

Pij = I3il3j + auej, 

(2.14) 

Pi68 = Pi (3/36 + Ps) + aneies, 

(2.15) 

P79 = 3/37 + P2 + Oh £79 

(2.16) 
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for i = 1, 3, 7 and j = 2,4, 5, 6 . The parameters f3i are independent of an- They are dehned 
in the appendices of Refs. [22, 23] and, hence, shall not be given here again. In Eqs. (2.13) 
to (2.16), we extend this parametrisation to theories beyond Horndeski gravity (q;h 7 ^ 0), 
introducing the parameters Cj dehned in Appendix B. Finally, we have also made use of the 
sound speed of the dark energy component [22, 23] 


c 


2 

S 


2(1 + Q!b)^ 
a 


1 + CKT 


1 + Oh 
1 + Ob 


^1 + am 



/l + aH\ (1 + Oh)^ Pm 

Vl + aB/ J a WIP' 

(2.17) 


where a = 6 ag + ax- 

A second held equation follows from using tt' determined from the time-time component 
Eq. (A.l) in the time-space component Eq. (A.2) to get vr. From the traceless space-space 
component Eq. (A.3), it then follows that 


+ kl{Bi^ + B2^) = Pm (^sAm + ^4^) , 
where the coefficients are given by 

(1 -|- ot) + Q!b) + ^^^"79 + 


Bi = 
B2 = 
Bs = 


54 = - 


19 + oiniaB - au)k‘jj 
+ («B - Qh)^ 

79 -F aH(aB - aH)k%' 

1 _ ^79 + QiHeB3 

79 -K aH(aB - an)k‘jj ’ 

1 ^79 + aHeB4 + aH(aB — aH)fc|/ 


252M2 


79 + aH(aB - OLB)k\ 


H 


(2.18) 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 
( 2 . 22 ) 


with 79 = 79 — and eBi dehned in Appendix B. The parameters 7 * are independent of 
an- They can be found in Refs. [22, 23] and are not given here again. 

Finally, in addition to Eqs. (2.8) and (2.18), we derive the energy-momentum conser¬ 
vation equations from the modihed Einstein equations in Appendix A, which will close the 
system of differential equations. Using the time derivatives of the time-space and trace of the 
space-space components to eliminate with tt' from the time-space component, we hnd the 
usual momentum conservation equation 


Urn + 


(2.23) 


Using all of the modihed Einstein equations and the scalar held equation together with the 
time derivatives of the time-time, time-space, and traceless space-space components to solve 
for A', A, vr and T and their time derivatives, as well as we derive the standard energy 
conservation equation 


Am = -knVra “ 3C', 


where we have used the comoving curvature (" = — Vra/kjj. 


(2.24) 
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3 Semi-dynamical linear cosmological perturbations 

As emphasised in Secs. 1 and 2, performing a quasistatic approximation of the modified Ein¬ 
stein and energy-momentum conservation equations, Eqs. (2.8), (2.18), (2.23), and (2.24), 
is not sufficient for conducting a more efficient but comprehensive test of gravity. The ap¬ 
proach is limited by its incomplete description at very large scales and cannot constrain the 
fully available EFT model space. An analysis beyond the simple quasistatic approximation is 
therefore needed, implying the consideration of velocities and time derivatives of the metric 
fluctuations. In order to account for the contribution of these terms in the field equations, we 
propose an extension of the quasistatic approximation that introduces corrections from 
and Vm/kn in Eqs. (2.8) and (2.18), determined at a pivot scale /c* of choice. More specif¬ 
ically, we solve the coupled system of modified Einstein and energy-momentum conservation 
equations at /c* and then use this result to extrapolate relations between the perturbations 
to other scales by performing the replacements 

^ (3.1) 

C' ^ (3.2) 

in Eqs. (2.8) and (2.18), where /$ = 4>'(A:*)/<h(A:*) and ^characterise the 

rates of evolution of the spatial metric potential and the comoving curvature at the pivot 
scale, respectively. Furthermore, from momentum conservation Eq. (2.23), we obtain 

^ ^ I 7 [(/<!> - k)^ - • (3-3) 

These replacements allow us to eliminate time derivatives of $ and the velocity field in the 
modified Einstein equations. 

Motivated by the energy-momentum conservation equations in Sec. 2 and the gauge 
transformation of the comoving curvature, which are not modified in unified dark energy, we 
estimate that relations between the fluctuations scale as 

(3.4) 

Hence, the field equations and the coefficients Ai and Bi suggest that the replacements pro¬ 
posed in Eqs. (3.1) through (3.3) will provide a good approximation for Horndeski scalar- 
tensor theories at small scales as the velocity field and the time derivatives of the metric 
potentials are subdominant and do not contribute in the limit of fc —>■ 00 . However, if we 

allow an 7 ^ 0, a dependency on at A; —)> 00 enters through Eq. (2.18) and a dependency 

on Vm/kn through Eqs. (2.8) and (2.18), suggesting that a quasistatic approximation with 
k = fc = ^ becomes invalid (cf. [39]). Hence, for beyond-Horndeski models nonzero contri¬ 
butions of /$ and should be considered at leading order. Note that here and throughout 
the paper, /c —)• 00 refers to a formal limit taken in linear theory. In concordance cosmol¬ 
ogy linear predictions are typically applicable to scales of A; < 0.1 /i“^Mpc. This limit can, 
however, be affected by modifications of gravity [40] and one should check against nonlinear 
corrections when using linear predictions to constrain the EFT model space with cosmological 
observations. 

Alternatively to the replacement in Eq. (3.2), one could define /a^hi = 

Eq. (2.24) to eliminate velocities in the field equations at small scales and derive a first-order 
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inhomogeneous nonlinear differential equation for /a from the energy-momentum conser¬ 
vation equations and Eqs. (2.8) and (2.18). This would provide a good approximation at 
small scales but not facilitate computations at large scales. For an approximation at super- 
Hubble scales (/c —)• 0), we could also use as defined in Ref. [ 8 ] to replace velocities, i.e., 
limfcj:^_,.o C = However, this comes at the cost of introducing an extra scale depen¬ 

dency in Eq. (3.3) that consequently adds further coefficients in the expansion of Eqs. (2.8) 
and (2.18) in powers of k and, hence, in the semi-dynamical description of effective unified 
dark energy modifications discussed in the following. 


3.1 Effective modifications 

Performing the replacements Eqs. (3.1) through (3.3) in the modified Einstein equations of 
the Fourier amplitudes of linear cosmological fluctuations, Eqs. (2.8) and (2.18), we obtain 


H \ H 

/4 + /I + ( “ ^AjpPmj /$ + — + A2k‘jj 




1 ~ -B4 —Pm I /$ -|- H4 —Pm/c + Blkjj + ( B^ — p^ -|- B 2 k‘fj ) 'k 




H 

^W" 


^3Pm^m(3.5) 

H3PmAm(3.6) 


Having cast the field equations in this form, we combine them and parametrise the gravita¬ 
tional properties of unified dark energy as an effective modification of the Poisson equation 
p and by the introduction of a gravitational slip 7 , 


^ 2H‘^k]j^ _ 1 P+2k]j + + p+efcg 

^^PmAm P—0 T P‘—2k‘^ P—4fc|^ -|- P—g/c^ 

^ ^ 7-ro + l+ 2 kH + l+Akjj 

'I' P-+2 + P+4k‘jj + p+ekjj 


(3.7) 

(3.8) 


respectively, where p± 2 n and 7 + 2 n {n = 0,1, 2, 3) are functions of time only and presented 
in Appendix B. In combination with energy-momentum conservation, Eqs. (2.23) and (2.24), 
and given /$ and /(^, these two relations completely describe the linear fluctuations of unified 
dark energy in the semi-dynamical approximation and constitute the core of the approach 
from which further relations can be derived. 

In the limit of A: —>■ 00 , Eqs. (3.7) and (3.8) simplify to 


_ p+6 _ 1 k‘to + “H {hpt,oc + /c/^Joo) 

p_6 k2m2 ^ + fcPlJ) 

^ J±± ^ _ 

pi + OH 


(3.9) 
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(3.10) 













where the time-dependent coefficients are given by 


/^oo 


f^oo 


(q;b — Oh)^ Q^b (1 + Q^t)/^? + CKH ^/ 3 ' 


VCBl — 2eB3CE 


H' 

+ 2 q;b (^7 - Cs )79 (■ ^ 






a(aB — Oi}iYf5T— — 2en(aQ;HeB3 + «b79) 

ri 

en(aaHeB3 + aB79) 


Pm 


2 i? 2 M 2 ^ 


Pm 


i72M2’ 


a{aB — O'uYcl— + anen {« [«b(1 + Q^t) + aneBi] + 2aB79} 


Pm 


2^2 M 2 ^ 


q;Q!b(ckb ~ ckh) 7 
—aen(aB — an) 


Pm 


2 iJ 2 M 2 ’ 
2 Pm 
2 iJ 2 M 2 ’ 
H' 


a{aB — an) + «Hen(aaHeB3 + aB79) 


Pm 




(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 


For Horndeski theories, where oh = 0, this further simplifies to (cf. [22, 23]) 


Poo — 

7 oo — 


1 2 [aB(l + CTt) ~ «M + C^t]^ + «(1 + CkT)Cs 

K2_/Vf2 aCg 

2aB [crB(l + ckt) ~ ckm + ckt] + CKCg 
2 [aB(l + «t) ~ «M + ot]^ + ct(l + a^)c^ 


(3.18) 

(3.19) 


As expected, while for Horndeski theories, poo and 700 are independent of /$ and /(^, this is 
no longer true if an Y Oj implying that the quasistatic approach /$ = /c; = 0 is generally not 
a good approximation for beyond-Horndeski theories. Hence, when an Y 0; even in the limit 
of /c —>• 00 , neither 4>' in Eq. (2.18) nor Vm/kn in Eqs. (2.8) and (2.18) should be neglected. 

It is also worth noting that a combined measurement of the cosmological background 
as well as poo and 7^0 through the large-scale structure only provides constraints on H and 
two ai functions. This generally leaves two Oj functions unconstrained in Horndeski gravity 
and an additional one in beyond-Horndeski models. In principle, one could obtain additional 
constraints on ctt and ctm from a measurement of the linear propagation speed and friction 
of gravitational waves [31, 34] but the required observations attributed to a late-time cosmo¬ 
logical effect are not yet available. Also note the independence of Eqs. (3.18) and (3.19) on 
ax- Hence, to fully constrain the available model space in the EFT formalism another two to 
three measurements are required, for instance of p and 7 at additional, larger scales, where, 
however, V^/kn and time derivatives of may not be neglected, and for which we envisage 
the use of the semi-dynamical expressions Eqs. (3.7) and (3.8). Importantly, a quasistatic 
consideration of these additional scales instead suffers from degeneracies between specific re¬ 
lations in the EFT coefficients and concordance cosmology and contains no information on 
M|(t) [28]. 


3.2 Growth of structure 

In the semi-dynamical approach the logarithmic growth rate of matter density fluctuations 
/a = can in principle directly be determined from energy conservation and the 

evaluation of /$ and at the pivot scale A;*. However, in the determination of /a, deviations 
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from the pivot-scale evolution become important since velocities contribute at first order to 
through energy conservation. More specifically, from Eq. (2.24) and the replacements in 
Eqs. (3.2) and (3.3), we find 


/a = -^7 


H 

IT 


/$ “ /c 4— 

7 


3/^1 

^ kjj \ 2m ’ 


(3.20) 


where the modification of the Poisson equation p, and the gravitational slip 7 are given in 
Eqs. (3.7) and (3.8), respectively. In the concordance model this reduces to the precise relation 
/a = 1 -|- /#. The accuracy of Eq. (3.20) can, however, suffer for models introducing a strong 
scale dependency in /$ and /^, rendering a pivot-scale evaluation insufficient. 

The logarithmic growth rate /a relates velocities to matter density fluctuations via 
energy conservation in Eq. (2.24). The proportionality between the two becomes exact when 
Q' — 7 > 0. To improve the accuracy in the recovery of this relation, instead of using energy 
conservation, we can adopt the parametrisation = —huVra [20]. This allows us to 

formulate an evolution equation for using momentum conservation in Eq. (2.23), 


T' + T‘^ + 


, H' 


K^Pm 


= 0 , 


(3.21) 


where we have used that and /a relate as 


= + (3.22) 

Note that and /a agree in the limit of A; —)> 00 and Eq. (3.21) reduces to the more familiar 
quasistatic differential equation for /a. Unlike that equation, however, Eq. (3.21) applies to 
all scales and is not limited to quasistatic scales only. 

The parametrisation can be motivated by its observational implications, describing the 
relation of redshift-space distortions to the underlying matter density fluctuations, as these 
measure velocities rather than growth of structure /a directly. One can then furthermore 
directly determine from the modifications [ 20 ] of the relativistic effects contributing to 
galaxy clustering [41-44], where a non-quasistatic treatment of perturbations can become 
essential. These deviations can be large for models like Dvali-Gabadadze-Porrati braneworld 
gravity [45], where 7 can be of order unity at the largest observable scales [ 8 , 46], or for linearly 
shielded modified gravity models [28], which recover ACDM on quasistatic scales. If measured, 
these effects may provide new insights to gravitational physics and dark energy near the 
Hubble scale. Finally, Eq. (3.21) may also serve to formulate a growth-index parametrisation 
of matter density fluctuations 7 [47, 48], or rather of the velocity-to-density ratio, through = 
f2m(«)^, where Hm(a) = K^/Om/(377^). This is a convenient parametrisation for observational 
tests of modified gravity and dark energy (e.g., [49]), and in this form can be made consistent 
for applications beyond quasistatic scales (cf. [50] ). For small modifications of the concordance 
model, Eq. (3.21) may also be linearised in the deviations to simplify computations [51]. 

3.3 Super-Hubble evolution 

While the description of the modified structure is simplified in the limit of small scales, 
where as discussed in Secs. 3.1 and 3.2 the approximation /$ = /(; = 0 becomes accurate for 
Horndeski theories to leading order, a similar approximation can be made at super-Hubble 
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scales. More specifically, in the limit of /c —)• 0, metric theories of gravity with energy- 
momentum conservation satisfy the adiabatic fluctuations of a flat universe with conservation 
of the comoving curvature = 0 [52, 53]. This has been shown to apply explicitly to the 
unified dark energy description used here in Ref. [23]. As a consequence, from momentum 
conservation it follows that 


n" n" / n' n"\ 

C" - —C' = d>" - T' - —- ( — - — ) T ^ 0 . 

^ H' \H H' 


(3.23) 


Hence, in this limit, we can set /( = 0 and determine the evolution of the metric potentials 
from Eq. (3.23) given a relation between $ and 'I' from the Einstein equations. In ACDM 
and dark energy models, where 7 = 1, the evolution of the Newtonian potential obtained 
from Eq. (3.23) applies on all scales and for the concordance case with ^ = 1 also determines 
the matter density fluctuations and velocities. Hence, for small deviations from ACDM, we 
expect that a pivot scale of /c* = 0 provides a good approximation for /$ and at A; > 0. 

We can also directly derive an evolution equation for /$ at /c —)• 0, adopting Eq. (3.4) 
to set /( — 7 - 0 and Am —?■ 0 in the modified Einstein equations, Eqs. (3.5) and (3.6), which is 
equivalent to Eq. (3.23) supplemented with a modified Einstein equation. This yields 


/4 + /I + 



/3i2 ^ 

h /3i 

l^ 

h 


2aj3i 

+ -F 

h 


+79):g§^ _ 

^ — (Q:Q!HeB4 + «K79) 
a(l aT)/3i687iHW 
2aPi^ - (aaHeB4 + 



= 0 . 


(3.24) 


The unified dark energy modifications /r and 7 are then determined by Eq. (3.24), = 0, 

and Eqs. (3.7) and (3.8). Note, however, since this implies /^ = 0 on all scales, for models 
where contributes significantly at fc > 0, as we will encounter in Sec. 4, setting the pivot 
scale at /c* = 0 may not accurately recover the unified dark energy fluctuations. 


3.4 Choice of pivot scale 

In the following, we discuss a few natural choices for the pivot scale A:*, for which the modified 
perturbed Einstein equations, Eqs. (2.8) and (2.18), and energy-momentum conservation 
equations, Eqs. (2.23) and (2.24), can then be solved to obtain /$ and /^. This determines the 
next-to-leading-order terms for the small-scale modifications introduced in Horndeski gravity 
and describes the leading-order modifications in beyond-Horndeski theories. Depending on 
application, one may decide for different choices of pivot scale or use multiple pivot scales to 
improve the accuracy over a wide range of scales. 

• Super-Hubble scales: As discussed in Sec. 3.3, the comoving curvature is conserved 
at A:* = 0. Hence, in this limit we can adopt /^ = 0 to simplify the modified Einstein 
equations and describe the evolution of /$ through Eq. (3.24). This approximation can, 
however, become inaccurate for models where contributes significantly at A; > 0. 

• Hubble scale: The Hubble scale describes the visible horizon and hence for observa¬ 
tional purposes it is sufficient to restrict approximations to A; > aH and determine /$ 
and from solving the modified Einstein and energy-momentum conservation equa¬ 
tions at A:* = Hq. 
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• Sub-Hubble scales: Given that close to the Hubble scale observations are limited by 
cosmic variance, one may also choose to solve the field equations at a sub-Hubble scale, 
e.g., fe* = 20Hq, to improve the accuracy at the observationally more interesting scales, 
where measurements are more significant. For models where the sound speed of dark 
energy is significantly smaller than the speed of light, the contributions of velocity and 
time derivatives of the metric potentials at scales below the visible horizon can become 
important [30]. In this case, a natural choice for the pivot scale would be the sound 
horizon /c* = Hq/cs- 

• Small scales: Finally, the quasistatic approximation with /$ = /^ = 0, applicable 
in Horndeski gravity in the limit of A: —)• oo, allows one to determine the effective 
modifications from ACDM analytically from Eqs. (3.18) and (3.19) without the necessity 
of solving the field equations at a pivot scale, but it can break down near the Hubble 
scale. Moreover, as pointed out in Sec. 3.1, one should be cautious when applying this 
approximation when uh 7 ^ 0 , in which case /$ and may give rise to leading-order 
contributions in A: to /i and 7 . Hence, a sub-Hubble pivot scale may be of interest in 
particular for beyond-Horndeski models, with evaluation of /$ and at A:* —>• 00 . 

4 Examples in dark energy and modified gravity 

We test the semi-dynamical approximation introduced in Sec. 3 against numerical solutions of 
the exact modified linearly perturbed Einstein and energy-momentum conservation equations 
for the different free functions ai (t) introduced in the unified dark energy formalism described 
in Sec. 2. For illustrative purpose, we set cosmological parameters to Planck 2015 values [54]. 
In specific, we use the matter density parameter = 0.308 and dimensionless Hubble 
constant h = 0.678, corresponding to Hq ~ 2.26 x 10-^ /i/Mpc. 

Before studying effective modifications from concordance cosmology in the semi-dynamical 
approach, we quickly check that it consistently reproduces ACDM phenomenology in the cor¬ 
responding limit. In ACDM we have 0 :^ = 0 Vi with = 1 and = 

Hi[n mn ^ + (1 — Dm)]. Using these relations in Eqs. (3.7) and (3.8) for jjL and 7 with the 
coefficients defined in Appendix B, we obtain ^ = 7 = 1 , irrespectively of the values for /$ 
and as expected. Hence, the metric potentials match, <I> = —4', and the standard Poisson 
equation relates the Newtonian potential to matter density fluctuations. We can find the 
evolution of the fluctuations from solving 


which follows from Eqs. (3.23) or (3.24). 

Next, we allow the contribution of a kinetic term and potential of the scalar field to the 
Einstein-Hilbert action in Sec. 4.1, introducing ax 7 ^ 0 and H 7 ^ PAacdMj which represents the 
linear perturbation theory of quintessence models. In Sec. 4.2, we study metric f{R) gravity 
as a representative example for models with a running Planck mass (om 7 ^ 0) and braiding 
(ob 7 ^ 0). Finally, in Sec. 4.3, we alternatively allow a deviation of the speed of gravitational 
waves from the speed of light (ctt 7 ^ 0) and the introduction of beyond-Horndeski terms 
(an 7 ^ 0) to the action. In all scenarios, we assume general relativistic initial conditions in 
the matter-dominated regime. 
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4.1 Kinetic contribution and scalar field potential 


We first study models with contribution of a scalar field kinetic energy (ok 7 ^ 0) and potential, 
which allows for a free expansion history H. This scenario is realised in quintessence models, 
which are represented by ax = 3[1 — nm(fl)][l + 'U^(o)] with w describing the equation of state 
of the dark energy component and om = ob = ckt = 0. Thereby, Eqs. (3.7) and (3.8) become 



/ 21 ^ Pm \ ^ 



_ _ 

sw) + (3 + f") 


Pm 

H2M2 


2H" 

H 


ft^aK + kjj 


(4.2) 

and 7 = 1, respectively, where we have replaced using (^" = (/^ + /(^/$)d>. Hence, while 
there is no gravitational slip in quintessence models, we observe a clustering of the dark energy 
on large scales when /( 7 ^ 0. In the limit /c —?• 00 , we recover the quasistatic results /r = 7 = 1, 
which thus do not account for the clustering effect at kn ~ 1- Note that as we have /^ = 0 
in the limit of A; —)• 0, setting the pivot scale at /c* = 0 to extrapolate to sub-Hubble scales 
cannot account for the clustering either, hence, suggesting the choice of a sub-Hubble /c*. The 
evolution of /$ at /c —>• 0, following from Eq. (3.23) or (3.24), is determined by Eq. (4.1) as 
in the concordance model. 


In Fig. 1, we show numerical results for the semi-dynamical approximation of the mod¬ 
ification of the Poisson equation /r(a, k) for quintessence models with dark energy equation 
of state w = —0.8 and w = —a. Thereby, we chose a pivot scale of /c* = 20i4o) which accu¬ 
rately recovers the modifications at the observationally interesting scales of /c > WHq that 
are determined from solving the exact perturbed modified Einstein and energy-momentum 
conservation equations [20]. Finally, note that for the models studied here, due to the weak 
scale dependence in /$ and the logarithmic growth rate /a can directly be computed from 
Eq. (3.20) to good accuracy, as we have checked numerically. 


4.2 Evolving Planck mass and braiding 


A running of the Planck mass is described by om 7 ^ 0 whereas braiding or mixing of the 
metric and scalar kinetic terms is described by a-Q 7 ^ 0. A scenario in which both parameters 
are nonzero is, for instance, realised in f{R) models, where the free nonlinear function f{R) 
is added to the Einstein-Hilbert action [55], and represented by om = 2aB = //?/(!+ /h) 

Ok = ckt = 0 with subscripts of R denoting derivatives with respect to the Ricci scalar. In 
f{R) gravity Cg = 1, which can be determined from Eq. (2.17) using the Friedmann equations 
that can be found in Refs. [18, 19, 22, 23, 28]. This simplifies expressions for fx and 7 to 


1 /i+2/c|^ + laMkjj 

k‘^M^ -I- fi-2k‘jj- + ^aukjj ’ 


7+0 + 

fi+2 + ^aMk%' 


fx{a,k) 


l{a,k) 


(4.3) 













Figure 1. Effective modification of the Poisson equation fi{a, k) determined from exact linear pertur¬ 
bation theory (solid lines) and the semi-dynamical approach (dashed lines) for a quintessence model 
with constant and time-dependent dark energy equation of state, w = —0.8 (left panel) and w = —a 
(right panel), respectively. A quasistatic approximation of the modified Einstein equations yields 
^ = 1 on all scales as in the concordance model. The semi-dynamical approach introduces corrections 
to this result that match the exact perturbations of the dark energy models and describe the clustering 
of the quintessence field at very large scales. We chose a pivot scale of fc* = 20iJo — 0.005 /i“^Mpc 
to cover the observationally interesting scales of A: > IOHq. 


where = 1 + /it and the coefficients of Eqs. (4.3) are given by 
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(,4.5) 

(4.6) 

(4.7) 


In the limit of A; —>■ oo, we recover p^o = (1 + /h)~^ 4/3 and 700 = 1/2. 

To test the performance of the semi-dynamical approximation against the exact per- 
tnrbations of f{R) gravity [56, 57], we set the pivot scale at A;* = 0 , where /^ = 0 and 
solve 
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Figure 2. Effective modification of the relation of the lensing potential ($ — 'I')/2 to matter density 
fluctuations E(a, k) = /r( 7 +l )/2 (left panel) and the gravitational slip 7 ( 0 , k) (right panel) determined 
from exact linear perturbation theory (solid lines), the semi-dynamical approach (dashed lines), and 
the quasistatic approximation (dotted lines) for f{R) gravity with |/_ro| = 0.01 {Bq ~ 0.05). The pivot 
scale is set at A:* = 0. Note that a conventional quasistatic approximation yields a scale-independent 
E = (1 -I- /fl)“^, corresponding to the fc —>■ 00 limit (fc « IOOTJq) in the semi-dynamical approach, and 
hence does not account for the scale dependence in S of the exact solution. 


which follows from Eq. (3.24). Note that Eq. (4.8) is equivalent to the super-Hubble evolution 
found in Eq. (43) of Ref. [57] with om = H'B / H, where the Compton wavelength parameter 
is defined as B = /h_r(1 + fR)~^R'H/H'. The resulting effective modifications are shown in 
Fig. 2 for a designer f{R) gravity model [57] with |/_ro| = 0.01 {Bq ~ 0.05) and an equation 
of state of the dark energy component of tc = — 1 , where subscripts of zero denote evaluation 
today. These parameter values are cosmologically not viable [58] and serve here only for 
illustration. To better illustrate deviations from the quasistatic approximation, we show 
S = |u( 7 -|-1)/2 rather than /r in Fig. 2, where S relates the lensing potential (4> — 'I')/2 to the 
matter density fluctuations using Eqs. (3.7) and (3.8). In the conventional quasistatic limit, 
neglecting velocities and time derivatives of and T in the field equations in Appendix A, one 
obtains a scale-independent S = (1 -|-/j?)~^ whereas the semi-dynamical approach introduces 
a scale dependency in S, however, with the effects being too small to be observed with 
lensing [59]. Finally, note that for f{R) models, due to the strong scale dependence in /$ 
and /^, determining the logarithmic growth rate /a directly from Eq. (3.20) does not yield a 
good approximation and one should integrate Eq. (3.21) to obtain more accurate results. 


4.3 Varying tensor speed and beyond-Horndeski terms 

Finally, we study scenarios in which alternatively ot 7 ^ 0 or op 7 ^ 0- We start with the 
tensor speed modification and set ot = Q;To[l~^m(fl)]- With a concordance model expansion 
history, we obtain 7 = kjj/{S ot f(^ + kj^) and /r = l + 3aTf(/k‘jj. Hence, for /( 1 deviations 

from the concordance model structure remain small at sub-Hubble scales. However, the 
expressions for 7 and // become more complicated when H deviates from its ACDM evolution. 
For illustrative purpose, we set re = —1/2 and oto = 0.1 and solve the modified Einstein 
equations in Appendix A for = 1 and Oj = 0 (i 7 ^ T). The super-Hubble evolution 
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Figure 3. Effective modification of the Poisson equation /x(a, k) (left panel) and gravitational slip 
7 ( 0 , k) (right panel) determined from exact linear perturbation theory (solid lines) and the semi- 
dynamical approach (dashed lines) for a modification of the tensor speed with ax 7 ^ 0. The pivot 
scale is set to fc* = Hq. 


follows from Eq. (3.24), which here becomes 
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However, since /^ = 0 at A: —)• 0, we set fe* at the Hubble scale. The corresponding results are 
shown in Fig. 3, where we find good agreement between the semi-dynamical approximation 
and the exact modifications. We also find that if choosing a pivot scale at sub-Hubble scales 
for the model adopted here, Eq. (3.20) provides a good approximation to the logarithmic 
growth rate /a since /$ and are only weakly dependent on scale within the visible horizon. 

Finally, we analyse the new phenomenological aspects that beyond-Horndeski theo¬ 
ries [37] introduce with ap 7 ^ 0. Thereby we adopt ap = aHo[l ~ f^m(o)] and require 
= 1 and aM = ok = ctT = 0. In order to have a well-defined general relativistic limit 
^ = 7 = 1 when ap —>• 0, we set ap = —an- For illustrative purpose, we set apo = 1/2 and 
also adopt w = —1/2 as in the case of ax 7 ^ 0, and solve the modified Einstein equations 
in Appendix A in the corresponding limit. The super-Hubble evolution of /$ reduces to 
Eq. (4.1). As discussed in Sec. 3, in case of nonzero an, and Vm/kn contribute to 7 and 
^ through /$ and at leading order in k. Hence, setting the pivot scale at fe* = 0, where 
/(^ = 0, does not account for the full correction. In Fig. 4, we compare the modifications // 
and 7 against the semi-dynamical approximation using the pivot scale k^ = IOOHq [k —>■ 00 ) 
instead, finding excellent agreement between the two for scales of A: > WHq. To illustrate the 
importance of nonzero contributions of /$ and /^, we also show the resulting modifications at 
k = lOOEAo derived in the quasistatic approximation of setting /$ = /^ = 0 in Eqs. (3.7) and 
(3.8). Moreover, we show results from a more conventional quasistatic approach of assuming 
domination of Am over V^/kn and neglecting time derivatives of the fields in the modified 
Einstein equations given in Appendix A in the limit of A; —>■ 00 , which yields 
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Figure 4. Effective modification of the Poisson equation ^(a, k) (left panel) and the gravitational slip 
7(0, k) (right panel) determined from exact linear perturbation theory (solid lines), the semi-dynamical 
approach (dashed lines), and quasistatic approximations for a modification with beyond-Horndeski 
term oh ^ 0 . The pivot scale is set to = IOOHq {k — >■ 00). The two quasistatic approximations 
are obtained by setting = 0 in the semi-dynamical (sd; dot-dashed curve) equations with 

k = IOOHq and by a conventional approach (cv; dotted curves) of neglecting V^/kn with respect to 
Ani and time derivatives of the fields in the modified Einstein equations in Appendix A in the limit of 
fc —>■ 00. The velocity field and the time derivative of the metric potential $ contribute to /r and 7 at 
leading order in k, which prevents the quasistatic approximations to reproduce the exact perturbation 
results. 


and 7 = 1 -|- ajj- In both approximations, there are 0(1) deviations from their counterparts 
computed from the exact perturbations, as expected (see Fig. 4). Note, however, that the 
quasistatic approximation within the semi-dynamical approach, setting f^ = f( = 0 , performs 
well at high redshifts, where laql ^ 1 , which is in agreement with a similar analysis performed 
in Refs. [39, 60]. Finally, we find that due to the weak scale dependence in /$ and /^, 
Eq. (3.20) provides a good approximation to the logarithmic density growth rate /a in the 
scenario discussed here. 


5 Discussion 

The development of a generalised description of the cosmological structure formed in modified 
gravity and dark energy models has been a very active field of research in recent years. One 
such formalism resulting from this exploration is the effective field theory of cosmic acceler¬ 
ation or unified dark energy [18, 19], encompassing Horndeski scalar-tensor theories [35] and 
beyond [37]. While the full system of differential equations in the EFT framework can be 
solved numerically, this can be computationally challenging as well as inefficient depending 
on application. Moreover, the model-specific solutions may only offer little insight into the 
generic behaviour of the wide range of modifications possible. Hence, a quasistatic approx¬ 
imation, neglecting time derivatives with respect to spatial derivatives, is often applied to 
simplify computations. This allows one to extract information about the relations between 
the fluctuations without the need of integrating the field and conservation equations. The ap¬ 
proximation can, however, break down near the Hubble scale or the sound horizon for models 
with subluminal sound speed of the dark energy, affecting tests of gravity and dark energy 
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employing observations at very large scales [20, 29, 30]. Furthermore, an analysis beyond the 
quasistatic regime becomes necessary to constrain the entire available EFT model space and 
break degeneracies with ACDM for particular subsets of this space [28], hence, demanding a 
more general, yet efficient approach. 

Here, we have developed a semi-dynamical treatment of the linear cosmological pertur¬ 
bations of unihed dark energy that accounts for the time derivatives of the metric potentials 
and large-scale velocity helds. Our approximation extrapolates the evolution of fluctuations 
at a pivot scale to arbitrary scales and facilitates the analysis of modihcations from the concor¬ 
dance model in its application to observational tests of gravity and dark energy. We adopted 
the unihed dark energy formalism of Ref. [22, 23] and derived the energy-momentum conser¬ 
vation equations and a set of combined modihed Einstein equations, eliminating the scalar 
held huctuations, for the fully available model space, including beyond-Horndeski models. 
Evaluating the velocity held and rate of growth of the spatial metric potential at the pivot 
scale, we determined the time and scale dependence of the effective modihcations from the 
concordance model in the form of a deviation in the Poisson equation and a gravitational 
slip and approximated the growth rate of matter density huctuations. At small scales, and 
for Horndeski theories, the resulting modihcations recover a quasistatic approximation while 
accounting for corrections at the next-to-leading order in /c, which become important near the 
Hubble scale. For models beyond Horndeski gravity, we found that the velocity held and time 
derivative of the spatial metric potential can generally not be neglected and contribute at 
hrst order even in the small-scale limit. Finally, we tested the semi-dynamical approximation 
against numerical solutions of the exact modihed Einstein and energy-momentum conser¬ 
vation equations for a range of dark energy and modihed gravity models, including kinetic 
contributions of the scalar held, running Planck mass, kinetic braiding or mixing, deviations 
between the speed of gravitational waves and the speed of light, and beyond-Horndeski terms. 
We found good agreement between the approximations and the exact results. 

Since we only require the integration of the held equations at a single scale, the semi- 
dynamical approximation can be used to signihcantly enhance the efficiency in numerically 
intensive observational tests of gravity and dark energy, or to perform analytical operations 
that help to gain more insights on the rich phenomenology of unihed dark energy. 
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A Modified perturbed Einstein equations 


We present here the modified perturbed Einstein equations in the unified dark energy for¬ 
malism of Refs. [ 22 , 23] that we have used in Sec. 3 to describe the effective modifications 
of the Poisson equation ^{a,k) and the gravitational slip 'y{a,k). Metric perturbations are 
described in the Newtonian gauge with the FLRW line element given in Eq. (2.2). However, 
contrary to Refs. [22, 23], we adopt the total matter gauge for the matter fluctuations. 

With these specifications, the time-time component of the modified perturbed Einstein 
equation, or Hamiltonian constraint, becomes 
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The traceless space-space component, or anisotropy constraint, is given by 
(aT + 1)^ + (®H + l)'k — uytHtt' -|- (q!m ~ Oi'-£)H'K = 0 
and the trace of the space-space component, or pressure equation, by 
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Finally, the scalar field equation becomes 
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B Coefficients 


Finally, we provide a list of the coefficients appearing in the combined perturbed modified 
Einstein equations, Eqs. (2.8) and (2.18), in Sec. 2 and contributing to the semi-dynamical 
approximation of fi{a,k) and 7 ( 0 ,/c) in Eqs. (3.7) and (3.8), respectively, in Sec. 3.1. Note 
that the parameters f3i and 7 ^, which are independent of an since defined for Horndeski 
models, are provided in Refs. [22, 23] and are not given here again. 


Combined modified Einstein equations The beyond-Horndeski terms (oh / 0) in the 
combined field equations and the modifications ^ and 7 are defined as follows: 


Cn =- (ob — an)" 

a 

_ Q 

ei - 

€2 = (/Ss + e4a)ei, 

2 

£3 = - 


Q^B _ _|_ Pm 

Ob ~ CKH V®B CKH 


64 = 


/34 — (1 + Ot) ( 1^2 — 


£2 

£l 


£5 = (“B — CKH — 1) 


£ 1 , 

H”' 


(B.l) 

(B.2) 

(B.3) 

(B.4) 

(B.5) 


h'Y 

n) ^ H 


+ 


Oh +1 
4a 


(12aB -OK- ban) ( ^ 2 ^ 2 ) 


1 

+ 2 


— — aB(aM — /?3 + 3) — 2(aM — ax — Q^h) — ci:h(q;m + 1) H—^ 
a an 




Pm 


+ ^ £5b + (Q^H + 1 ) 


11 , 

- -I— (6aB — Ok — 3aH|aB + Ij) 
2 a 


£6 = /3i £7 — — [(ob — Q;h)£3 + /32 — /^s] f — 2 <[ (1 -I- ob) 


772 m 2 

+ £50 (B- 6 ) 


Pm ] H' 

H 

a]i\H' 

2 - ] ~F7 ~ “M + CUT 

CXB J H 


a 


—ohob ( ^ I + 


a 


H 


a 


B 


a 


H 


2aB — Oh Pm \ £1 


£7 = 


2 aB 


a 


Ob 

Pm 


2i72M2 J 


a 


C^B CUR J Ob 27727^^2 


£168 = 3 


2agaK — obQ^k ^ , n 

-Ob ( £4a ~ ClM H- 77 + 2 


a 


H 


om ~ (St 




o 12 / / H' H” 

-3aB ( (S^ — OMfST + (Sm — (Sx — -^ 


£79 — 


Gob 


a 


a' 


a 


2 1 ^ ^ - £ 4 a 


£bi-(1 + «h) + 

1 (H' 

£B 3 = 7 TT + 


(SB (SH 
H' 

I 

Pm 


1 


Pm 


(SB — (SH 
— (Sb(1 + (Sm), 


a 


a 


2 ( _ rH 


(SB (SH 


H' 

-I- OM - (Sx - ^ 


2\H 2i72M2 


H' 

£B4 = —3aB^, 


(B.7) 

(B. 8 ) 

(B.9) 

(B.IO) 

(B.ll) 

(B.12) 

(B.13) 


- 20 - 



































_ «B H' 1 (H' \ 


(B.14) 


esa = GaBOH 


«M-(«H + 1) — CKT — Oh + 

as 


3(aM + an) — 2aB { 1- 

an 


/ \ / 

— CiR + 2) H-20^ + 2aMOH + 1 

an 


+ Ok 
— (an + 1 )ok) 


a 


av 


£56 = 2(aB + 1)— — aB(3aM + 2aT — /^s + 5) + (an + 1) 
an 


. Qi'c 


3(q!m + 1) ~ 2^^ — /^s 
Ob 


—2(aT + 1), 

e5c = 2 [Q:M(aH + 1) — ax + an] 


(B.15) 


(B.16) 


CXt. 


as 


+(aB — an) 


2 ( OM - CtT - ^ j - (GaBttH + «k) 


a' 


H 

OH 


+aB 


Ol 

Q!t — «T + (®M + CKT + 3)(aM ~ /^s) + /^S ~ 2 (q;]V[ — 2aT — 1)^^ ~ Ci'yi + 2 


an 


— < 2(q;m ~ ax)^^ + 2a3(ax + 1) + an [aM(aM ~ /^s + 2) — — /Js + l] 


an 


—otyi + 2aH(“T + 1) + ax + ax(/33 “ 1) + aM(aM — ax — /^a + 1)| • (B-17) 

Semi-dynamical coefficients We provide here the coefficients of the semi-dynamical mod¬ 
ifications p{a,k) and 'y{a,k), Eqs. (3.7) and (3.8), respectively, in Sec. 3.1. The numerator of 
p{a, k) and denominator of 7 ( 0 , k) are determined by 
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where we have defined /r+g = pF, + an (yf^pt>oo +/c/^^ 00 ) ™ ^9®- (3-11) through (3.13). 
The denominator of p{a, k) is given by 
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where we have defined ^_6 = /x^o + (^f^P^oo +/c/^Ccxd) ™ (3-14) through (3.16). 

Finally, the numerator of 'y{a,k) in Eq. (3.8) is given by 

7+0 = OiPiQPi— + 2(/3i 68 + ^i6)(a«HeB3 + «b79) 2H^m‘^ ’ (B.23) 
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where we have defined 7+4 = 7 ^^ in Eq. (3.17). 
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